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Abstract. The development of surveys which will be able to cover a large region 
of the sky several times per year will allow the massive detection of transient events 
taking place in timescales of years. In addition, the projected full digitalization of the 
Harvard plate collection will open a new window to identify slow transients taking 
place in timescales of centuries. In particular, these projects will allow the detection 
of stars undergoing slow eruptions as those expected during late helium flashes in the 
post-AGB evolution. In order to identify those transients which correspond with late 
helium flashes the development of synthetic light curves of those events is mandatory. 
In this connection we present preliminary results of a project aimed at computing grids 
of theoretical light curves of born again stars. 



1. Introduction 

About 20% of the stars leaving the Asymptotic Giant Branch (AGE) are supposed to un- 
dergo a final thermal pulse during their post-AGB evolution (Iben.il984). In those cas es. 



the star is predicted to be temporarily re born as a ye llow giant ( Schoenbemejl 1979h in 



the so-called Born Again AGB scenario (!lbenlll984h . The transition from the prewhite 



dwarf to the giant configuration is exp ected to be very rapid; bein g of a few years in the 



very late thermal pulse (VLTP) case dlben & MacDonald 1995) and of the order of a 
century in the case of a late thermal pulse (LTP) (" Schoenbemer 1979 Y Due to the short 



duration of these events, although ~ 10% of all post-AGB stars are expected to undergo 
a VLTP, the observation of those events is very rare. Indeed, only two objects (V605 
Aql and V4334 Sgrfl have been identified as stars undergoing a VLTP ( Duerbeck et al.l 
2000L 12002). while a third has been identified with the LTP flavor of the scenario (FG 



Sge, see iJefFery & Schonberner 2006 and references therein). Individual born again 
stars are very valuable to test our understanding of the s-processes during the thermal 



'In addition two other erupt ions have been suggested to be VLTP events; CK Vul dEvans et al]|2002l) and 
more recently NSV 1 1749 I IMiller Bertolami et alj20llb . but neither of these have been confirmed. 

1 



2 



Marcelo M. Miller Bertolami & Rene D. Rohrmann 



pulses on the AGB ( Asplund et all 1 19991 : 1 Jeff ery & Schonbemeill2006l) and also to test 
and validate iJie sim ulations of convective reactive burning in the interior of the stars 
dHerwig et al.ll201lh . In addition, these events are key to un derstand the formation o f 
carbon-rich H-deficient stars and H-deficient white dwarfs ( Werner & Herwig]|2006h . 
Consequently, the identification of more stars undergoing late helium shell flashes will 
boost our understanding of the late phases of low and intermediate mass stars and their 
nucleosynthesis. In this connection, a rough estimate sug gests that the birth rate of 
planetary nebulae in our galaxy is of about one every year dZijlstral [20021) . If 10% of 
their central stars become VLTP giants, then we should expect such events in our galaxy 
to take place at a rate of about one per decade. Due to their high intrinsic brightnesses 
(Mv=-2...-4) these objects can be easily detected at large distances within our galaxy. In 
this context, the development of surveys which will be able to cover a large region of the 
sky several times p er year, like the Pa noramic Survey Telescope & Rapid Response Sys- 
tem (PanSTARR S,|lonrY^eLay|2M^ and the Large Synoptic Survey Telescope (LSST, 
IWalkowic3l2012 ). will allow the systematic detection of born again events. In addition, 
the projected full digitalizatio n of the Harvard plat e collection (Digital Access to a Sky 
Century at Harvard, DASCH. lGrindlav et al.ll2009h will open a new window to identify 
slower born again events taking place in timescales of centuries like FG Sge. Due to 
the lack of enough observational counterparts of bonafide born again stars the identifi- 
cations of such events has to rely partially in our present theoretical knowledge of their 
evolution. We have recently started a project to compute, and provide, theoretical light 
curves for born again events that will help in the identification of those transients. In 
what follows we present some preUminary results and discuss future developments. 



2. Numerical modeling 



In order to compute synthetic light curves appropriate for Bom Again AGB stars proper 
stellar evolution models and atmospheres are needed. As a first step we computed the 
evolution of stellar models during fast born again events (low mass VLTP events). Dur- 
ing VLTP events the hydrogen envelope remaining on the hot white dwa rf is convec- 

tively mixed into the interior of the star and violently burned — H-flash, see lMiller Bertolami et al. 
(12006). Stellar evolution sequences were computed with LPCODE stellar evolution code 
which is perfectly s uited to compute the mixing & burning during the violent H-flash 
(lAlthaus et al. 2005). Initial models for a wid e range of hot- white dwarf masses were 
taken from lMiller Bertolami & AlthausI (l2007b . which is the only available grid of de- 
tailed pre-bom again models to date. The masses of the models studied in the present 
work are 0.530Mo, 0.542Mo, 0.561Mo, 0.564Mo, 0.584Mo and 0.609Mo. As men- 
tioned above, in order to compute appropriate light curves for different filters, a proper 
assesment of the emergent spectral energy distribution of the models is necessary. Syn- 
thetic magnitudes for the stellar evolution models are computed with a modified version 
of the stellar atmosphere code described in Rohrmann et al. (2011) which includes a 
simplified treatment of the effects of C, N and O in the structure of the atmosphere and 
radiative transfer equations. In brief, the modelli ng of the stellar atm osphere includes 
several ions and molecules of H and He (see IRohrmann et al1l20I ll for a detailed de- 
scription) and the neutral and first three ionization stages of C, N and O. The chemical 
balance of these species is deterrn ined by the occupation probability formalism for H 
and He dHummer & Mihalaslll988 ) and by the ideal partition functions for the C, N and 
O ions. Opacity sources for H and He ions and molecules are included in a very detailed 
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and exhaustive way dRohrmann et all 1201 1 ) while C, N and O ion opacity sources are 
included within the hydrogenic approximation and no molecules have been included so 
far. 



3. Preliminary results and future work 

As expected, eruption light curves are dependent not only on the timescales set by 
the evolution of the stellar interior, but also on the specific element abundances of the 
photosphere of the models. In Fig. [T] (top right panel) we show the predicted light 
curves obtained from the same stellar evolution sequence when different surface abun- 
dances are adopted in the computation of the synthetic magnitudes. Preliminary syn- 
thetic light curves were computed for the UBVRJ (Bessell 1990) and the Sloan Digital 
Sky Survey (ugriz; Fukugit a et al.ll 19961) photometric systems. In Fig. [T]we show the 



computed light curves for all our stellar evolution sequences in the SDSS phot omet- 



ric sy stem, when abundances similar to those observed in V4334 Sgr (Asplund et al. 



are assumed. Note that the S DSS photometric system is similar to the one 



adopted by the PanSTARRSl survey d Tonry et"al]|20I2 ). In order to assess the im- 



portance of the detailed treatment of the atmosphere and also the importance of the 
different atmospheric abundances, two test cases have been computed by assuming dif- 
ferent atmospheric compositions for the 0.56 IMq sequence: a solar-like atmosphere 
iXH/XnJXc/Xt^/Xo/ = 0.74/0.25/2.4 x lO^VV x 10^4/5.8 x 10~^ brown lines) and 
a V4334 Sgr-like atmospher e (Xi^/Xij^JXc/Xn/ Xc,/ - 7.4 x 10^4/0.74/0.22/8.2 x 



10 ^/2.9 X 10 ^, black lines. lAsplund et ahll 19991) . Clearly, a correct assessment of the 



radiative transfer in the atmosphere is necessary to compute born again light curves. 

Several improvements in theoretical lightcurves, as well as a proper assessment 
of the involved uncertainties, are still necessary. In particular we plan to improve the 
stellar atmospheric modelling by including a better treatment of the C, N and O ions and 
the inclusion of related molecules. Also a calibration of the mixing efficiency in stellar 
evolution sequences during the violent H-flash is necessary and will be performed. In 
addition we plan to expand the grid of masses to include higher masses, with slower 
VLTP events, and to compute lightcurves for LTP episodes. 
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Figure 1. Synthetic light curves computed for the SDSS photometric system 
under the assumption of V4334 Sgr-Uke compositions for all our computed stel- 
lar evolution sequences. The effect of different compositions is shown in the top 
right panel for the 0.561Mo where the light curves for a solar-like atmosphere 
are also shown (grey curves). Note that the much slower eruption light curve of 
the O.6O9M0 sequence is due to the fact that massive remnants do not harbour 
enough hydrogen for the H-flash t o drive a fast expansion as already explained by 
iMiUer Bertolami & AlthausI (l2007h . 
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